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ABSTRACT 

Gradients  in  the  index  of  refraction  caused  by  temperature  fluctuations  create  optical 
turbulence.  This  optical  turbulence  greatly  degrades  the  performance  of  imaging  and  laser 
systems  operating  in  the  atmosphere.  There  are  several  methods  of  measuring  optical 
turbulence  in  the  atmosphere,  but  few  give  a  high  resolution  profile  of  the  turbulence  at  all 
levels  displaying  the  "layers"  of  turbulence  that  exist. 

This  thesis  is  an  investigation  of  the  use  of  a  high  speed  temperature  probe  coupled  with 
a  VIZ  9000  meteorological  system  to  measure  the  optical  turbulence  in  the  atmosphere  from 
a  ground  base  to  approximately  25  km.  The  processing  of  the  data  is  done  by  a  program 
developed  as  part  ol  this  research.  The  outputs  of  the  program  are  r0,  0O,  and  Cn:  as  functions 
of  altitude. 
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I.    INTRODUCTION 

Since  the  Space  Defense  Initiative  sparked  a  push  to  use  laser  weapons  in 
space,  many  scientists  have  devoted  their  attention  to  the  effects  of  propagating 
a  laser  beam  (high  power)  through  the  atmosphere.  The  problems  of 
propagating  a  laser  beam  in  the  atmosphere  can  be  broken  into  four  basic 
categories  of  detrimental  effects  [Ref.  1].  They  are:  absorption  and  scattering  by 
atmospheric  aerosols,  stimulated  Raman  scattering,  thermal  blooming,  and 
optical  turbulence.  It  is  important  to  note  that  all  of  these  effects  depend  not  only 
on  the  state  of  the  atmosphere  but  also  on  the  wavelength  of  the  propagating 
beam  and  the  mode  of  laser  operation  (CW  or  pulsed).  Significant  absorption 
of  incident  radiation  induces  thermal  blooming  when  the  absorbed  energy  alters 
the  atmospheric  index  of  refraction.  At  certain  wavelengths,  such  as  around 
3.8  [im,  atmospheric  "windows"  exist  where  absorption  is  lower.  Raman 
scattering  is  significant  only  when  the  peak  power  is  sufficient  to  induce  a 
nonlinear  index  of  refraction  in  atmospheric  nitrogen  and  oxygen,  inducing  higher 
order  harmonics  of  the  propagating  wavelength  [Ref.  1]. 

At  shorter  wavelengths  atmospheric  turbulence  and  scattering  tend  to 
dominate  the  beam  propagation  problem.  Turbulence  scrambles  the  phase  of 
a  laser  beam  across  the  wavefront,  thereby  spreading  the  energy  over  a  larger 
solid  angle  (^radians  to  milliradians)  than  would  be  expected  for  a  diffraction- 


limited  optical  system.  Scattering  from  aerosols  spreads  the  energy  over  very 
large  solid  angles,  on  the  order  of  one  or  more  radians,  which  severely  degrades 
collimation  and  reduces  the  beam  energy  that  can  be  delivered  to  a  small  target. 
When  light  propagates  through  a  medium,  the  local  phase  of  the  wave 
depends  on  the  product  of  the  index  of  refraction  of  the  medium  (n)  and  the  path 
length  between  the  source  and  the  point  in  question.  It  is  known  that  the  index 
of  refraction  of  our  atmosphere  at  a  given  wavelength  is  not  constant  because 
atmospheric  density  is  not  constant.  Some  contributing  factors  are  the  decrease 
of  density  with  altitude,  differential  heating  by  the  sun,  turbulent  wind  flow,  and 
convection  that  causes  localized  inhomogeneities  in  the  index  of  refraction,  called 
optical  turbulence  (small  variations  in  n).    The  index  of  refraction  in  air  is 


-6      P 


n  =   l+79xl(r6x—  (1) 

T 


where  P  is  atmospheric  pressure  (millibars),  T  is  temperature  (kelvin),  and  the 
humidity  term  has  been  ignored  since  its  contribution  is  so  relatively  small 
[Ref.  2],  Consequently,  light  that  originates  as  a  plane  wave  becomes  a  complex 
distorted  wave  after  propagating  a  finite  distance,  as  shown  in  figure  1-1  [Ref.  3]. 
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Figure  1  -1 :  Effect  of  Atmosphere  on 
a  Plane  Wave 


There  are  several  systems  designed  to  measure  optical  turbulence. 
However,  there  is  a  need  for  a  system  that  will  produce  a  high  resolution, 
accurate  vertical  profile  of  the  turbulence  verses  altitude. 

This  thesis  continues  the  development  of  a  single  point  temperature 
thermosonde  system  to  measure  the  vertical  profile  of  the  temperature  structure 


parameter. 


II.    BACKGROUND 

A.      THEORY 

From  the  arguments  in  the  previous  chapter,  it  is  apparent  that  analysis  of 
a  laser  system  that  propagates  a  beam  from  ground  to  space,  or  vice  versa, 
requires  an  accurate  understanding  of  atmospheric  turbulence  effects.  Although 
these  effects  are  random  and  complicated,  statistical  quantities  can  describe 
many  of  the  relevant  phenomena.  A  successful  method  that  has  been  developed 
for  such  analysis  uses  structure  functions  to  quantify  the  statistics  of  atmospheric 
turbulence  (rather  than  autocovariance  or  autocorrelation  functions).  Wavefront 
phase  distortion  depends  on  localized  differential  variations  in  n,  rather  than  upon 
those  properties  that  are  common  to  all  portions  of  the  wavefront.  Therefore, 
structure  parameters  are  inherently  differential  in  form.  A  structure  tensor  as 
defined  by  Tatarski  is  [Ref.  2]: 


Dij(r1,r2)=<[(f(r1)-f(r1))-(f(r2)-f(r2)  )]2>  <2) 


where  the  angled  brackets  represent  an  ensemble  average. 


By  making  three  assumptions  [Ref.  4]: 

1-  Isotropy,  so  that  only  the  magnitude  of  r  is  important,  not  the  spatial 
orientation  of  the  measurement 

2-  Homogeneity,  which  implies  that  a  change  in  position  does  not  alter  the 
statistical  properties  of  the  turbulence 

3-  Incompressible  turbulence 
then: 

Df  =  <  [fir^-f  (r2)]2>     ,  (3) 

and 


i  =    |  r1-r2  (4) 


Furthermore,  by  the  Kolmogrov  theory  of  turbulence,  as  long  as  r  is  within  the 
"inertial  subrange  of  turbulence"  [Ref.  5],  then 

Df  =  Cf2r2/2  (5) 

where  Cf2  is  a  structure  parameter.  This  inertial  subrange  of  turbulence  is  the 
range  from  the  small  scale  l0  of  turbulence  to  the  outer  scale  L0.  The  small  scale 
corresponds  to  the  eddy  size  where  viscous  dissipation  increasingly  consumes 
the  energy  in  the  eddy.  The  outer  scale  corresponds  to  the  largest  scale  size  for 
which  eddies  can  be  considered  still  isotropic  (see  figure  2-1) [Ref.  4]. 
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Figure  2-1 :  Hierarchy  of  Eddies  for 
Energy  Cascade 


For  temperature  fluctuations  that  follow  a  Kolmogrov  spectrum  the  temperature 
structure  parameter  CT2  is  [Ref.  2] 


2  <[T(I2)-T(l1))2> 


Ct  = 


^2-^1 


2/3 


(6) 


where  T  is  temperature  (kelvin)  and  r  is  distance  (meters). 

Similarly  the  index  of  refraction  structure  parameter  is  [Ref.  2] 


C  z  = 


<  [nir^-nir,)}^ 


^2-ri>" 


2/3 


(7) 


where  n  is  the  index  of  refraction  and  r  is  distance  (meters).  Since  atmospheric 
turbulence  pressure  fluctuations  are  small,  the  pressure  term  of  the  partial 
derivative  of  equation  (1)  is  negligible,  and  Cn2  and  CT2  are  related  by  the 
temperature  partial  derivative  [Ref.  6]: 


C^   =    (79xl(T6—  )2Cr 
T2 


(8) 


There  are  two  statistical  parameters  that  characterize  the  effect  of  optical 
turbulence.     The  first  is  r0,  which  is  a  coherence  length.     It  is  a  measure  of 


transverse  spatial  coherence  of  wave's  electric  field  [Ref.  7].  In  the  absence  of 
turbulence,  r0  is  comparatively  large  (see  figure  2-2).  The  importance  of  rG  is  that 
it  describes  the  spatial  resolution  of  an  optical  system.  It  is  often  called  the 
"seeing  parameter"  or  "Fried's  parameter",  and  it  determines  the  effective 
diffraction  limit  of  an  optical  system.  A  large  r0  is  desirable  if  resolution  is  to  be 
determined  by  the  characteristics  of  the  radiating  aperture,  rather  than  by 
turbulence. 

In  the  absence  of  turbulence,  the  resolution  of  a  diffraction-limited  system 
is 


eD=  -^    ,  o) 


whe.'e  6D  is  the  angular  divergence  of  the  beam  (in  radians),  A  is  the  wavelength, 
D  is  the  aperture  size,  and  a  is  a  geometric  factor  of  order  unity.  Beam  spread 
decreases  with  increasing  aperture  size  D,  until  D  reaches  r0.  Further  increase 
of  D  beyond  r0  does  not  result  in  decreased  beam  divergence;  the  system's 
angular  resolution  stays  constant. 
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Figure  2-2:  Symbolic  Projection  of  r 


Typical  values  for  r0  are  5  cm  (day)  to  10  cm  (night)  for  A  =  500  nm  [Ref.  8].  For 
a  plane  wave,  r0  depends  upon  the  path  integral  of  the  index  of  refraction 
structure  parameter,  Cn2.  From  Ref.  6 


■c  =  2.1  [1.46k2fc*(z)  dz] 


■3/5 


(10) 


where  k  is  the  wave  number.  Equation  10  shows  that  r0  is  proportional  to  km. 
Another  related  parameter  is  the  isoplanatic  angle,  90.  Whereas  r0  is  the 
spatial  coherence  length,  60  is  the  angular  coherence  length.  It  defines  the 
"cone"  in  the  atmosphere  where  turbulence  is  effectively  the  same  for  different 
angular  paths.  As  with  r0,  an  optimized  optical  system  should  be  characterized 
by  a  large  60.  60  can  be  found  by  two  methods.  First  it  can  be  calculated  by 
knowing  Cn2  along  an  optical  path  [Ref.  9]: 


0O  =    [2.9  0Sk2fc2n(z)  z5/3dz] 


■3/5 


(ID 
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where  z  is  a  path  length  in  meters,  or  it  can  be  measured  by  scintillation  of 
starlight  caused  by  random  variations  of  local  irradiance  induced  by  turbulence 
[Ref.  4].  By  use  of  equations  5  and  6,  a  balloon  system  with  a  high  speed 
temperature  probe  can  measure  a  profile  of  Cn2  from  the  ground  to  approximately 
20  km  above  the  earth.  This  profile  can  then  be  used  to  calculate  90.  Comparing 
optically  measured  60  values  to  those  computed  from  a  balloon  system  will  verify 
the  accuracy  of  these  measurements  [Ref.  10].  60  also  has  a  a6'5  wavelength 
dependence.    Typical  values  are  3-10  ^radians  at  X  =  500  nm  [Ref.  8]. 

Both  r0  and  60  are  measurable  parameters  that  determine  the  effects  of 
atmospheric  turbulence  that  degrade  the  collimation  and  spatial  coherence  of  a 
propagating  electromagnetic  wave.  However,  r0  and  60  are  independent  of  each 
other.  As  equation  10  shows,  r0  is  equally  weighted  along  the  beam  path.  Q0  is 
an  indicator  of  higher  altitude  turbulence,  however,  as  the  z5/3  weighting  in 
equation  11  suggests. 


B.       SYSTEMS  REQUIREMENTS 

1.       Processing  System 

In  1989  the  Naval  Postgraduate  School  began  development  of  a  single 
point  probe  thermosonde  system  to  obtain  a  CT2  verses  altitude  profile.  The  main 
processing  component  of  the  system  was  a  VIZ  9000  meteorological  processing 


system  that  utilizes  VIZ  9000  digital  Microsonde  II  thermosondes  attached  to 
weather  balloons.  The  earlier  version  sonde  (1543-364  Revx6)  has  8  channels 
to  collect  data,  but  only  4  of  these  channels  are  used  by  the  VIZ  system  for 
meteorological  data  collection.  The  remaining  unused  channels  are  used  by  the 
NPS  group  to  measure  the  temperature  of  the  atmosphere  with  a  high  speed 
thermistor  probe;  however,  only  3  of  these  4  channels  were  available  for 
measuring  temperature  due  to  manufacturer  design.  Newer  version  sondes 
(1 543-1 68B)  now  have  all  four  spare  channels  available  to  measure  temperature 
with  an  added  probe. 

The  earlier  NPS  single  point  temperature  probe  systems  used  a  0.5  x 
0.5  mm  square  Thermometries  thermistor  flake.  Since  this  thermistor  proved  to 
be  too  slow,  this  thesis  concentrated  on  a  faster  probe  [Ref.  10].  A  5  mil  VECO 
thermistor  probe  (page  10)  was  connected  to  the  thermosonde  PC  board  as 
shown  in  figure  2-3. 


CD4053  IC  CHIP 


16      15      14      13      12      11      10      9 
TO         •••••• 


PROBE 


L]       J      1       J      J    »  +  LEAD  TO 

PROBE 

SPARE  4 
SPARE  1 
SPARE  3 
SPARE  2 
5  VOLT  SUP 
Figure  2-3:  Connecting  Probe  To  Sonde 
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The  VIZ  9000  samples  8  channels  every  1.23  ±  0.01  second.  These  include  a 
calibration  signal,  pressure,  temperature,  relative  humidity  and  four  spare 
channels  (1  through  4).  With  the  exception  of  the  pressure  channel,  all  of  the 
measurements  are  resistances.  The  sampling  time  of  each  spare  channel  varies 
from  32  to  53  ms.  The  variable  time  is  a  consequence  of  an  algorithm  built  into 
the  MICROSONDE  II  sondes  that  compensates  for  the  104  to  107  range  of 
resistance  values  normally  encountered.  The  resistances  measured  by  the 
system  on  the  unused  channels  are  stored  into  a  computer  file  at  the  end  of  the 
flight  and  then  are  converted  into  temperatures. 

System  noise  measurements  of  the  older  version  VIZ  microsonde  II 
were  done  by  NPS  November  1990  [Ref.  10].  Table  1  is  a  summary  of  the 
measured  system  noise  for  the  new  version  sondes.  Fixed  resistors  were 
measured  by  the  sonde  giving  total  system  noise  expressed  in  terms  of  the 
standard  deviation  of  the  resistances  and  equivalent  RMS  temperatures. 


TABLE  1.  VIZ  MICROSONDE  II  SYSTEM  NOISE 

Resistor 

a  (ch  1)  (Q) 

a  (ch  2)  (Q) 

a  (ch  3)  (Q) 

a  (ch  4)  (Q) 

10  kQ 

3.9 
=  0.005  K 

3.7 
=  0.005  K 

3.5 
=  0.005  K 

3.6 
=  0.005  K 

43  kQ 

12.0 
=0.008  K 

11.8 
=0.008  K 

11.1 
=  0.008  K 

10.9 
=  0.008  K 

100  kQ 

21.8 
=  0.007  K 

20.7 
=0.007  K 

21.1 
=0.007  K 

22.4 
=  0.007  K 

1  M 

345 
=  0.012  K 

351 
=  0.012  K 

328 
=0.012  K 

333 
=  0.012  K 
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43  kQ 

12.0 
=0.008  K 

11.8 

=0.008  K 

11.1 
=  0.008  K 

10.9 
=  0.008  K 

100  kQ 

21.8 
=0.007  K 

20.7 
=  0.007  K 

21.1 
=0.007  K 

22.4 
=  0.007  K 

1  M 

345 
=  0.012  K 

351 
=0.012  K 

328 
=  0.012  K 

333 
=  0.012  K 

2.       Thermistor  Probe 

After  trying  many  different  probes,  the  NPS  group  has  had  the  most 
success  with  a  5  mil  VECO  thermistor  probe.  Figure  2-4  is  a  diagram  of  the 
construction  of  the  probe.  The  response  time  of  the  probe  was  found  to  be 


-  =  A  +  B/V 


(12) 


where  t  is  the  response  time  (sec).  A  and  B  are  constants  (A  =  0.00758  msec"1 
and  B  =  0.0077  msec"1  -sec0 5/m05),  and  v  is  the  wind  velocity  (m/s)  across  the 
probe. 
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Figure  2-4:  NPS  Single  Point  Probe 
Construction 
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Figure  2-5  shows  1  /t  verses  v° 5  for  wind  velocities  of  0,  1 ,  2  and  4  m/s  measured 
in  the  lab  with  a  chopped  Halogen  lamp  illuminator  as  a  heat  source.  The 
response  time  for  the  probe  for  an  average  ascent  rate  of  4  m/s  was  45  ±  5 
msec. 

To  calibrate  the  probe,  a  method  similar  to  the  Steinhart-Hart  equation 
[Ref.  11] 


-  =  a+Z?(ln  — )  +  c(ln  —  )2+d(ln  —  )3  (13) 

T  r0  r0  r0 


was  found  to  work  best  by  NPS  in  June   1990  where  a,   b,   c,   and   d  are 
coefficients  and  r0  room  temperature  resistance  [Ref.  12].    The  method  uses  an 
in  flight  calibration  of  the  probe  with  the  VIZ  temperature  sensor  as  a  reference. 
A  least  squares  algorithm  computes  the  coefficients  that  are  used  to  convert  the 
thermistor  resistances  to  temperatures. 
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Figure  2-5:  Probe  Response  Time 

The  structure  function  used  to  calculate  CT2  in  equation  6  is  sensitive 
to  the  response  time  of  the  probe  and  the  32  -  53  msec  integrating  period  of  the 
A-D  converter  that  is  used  in  the  VIZ  system.  Therefore  the  structure  function 
over  a  limited  scale  length  becomes  [Ref.  11]: 


sm(kxJ2) 


D(rtxvxJ  =4(sin2(— )/r5/3[ — ^-^]2[ ]dk 

{         2;  (kxJ2)    ,ll+(wt2)2J 


(14) 


where  a  is  2tt/U,  the  limiting  lower  frequency  (outer  scale  L0),  b  is  2tt/I0,  the 
limiting  upper  frequency  (inner  scale  l0),  k  is  2n/k,  the  wave  number,  t,  is  the  A-D 
integrating  time,  a>  is  kV  where  V  is  the  wind  velocity  (thermosonde  ascent  rate), 
and  t2  is  the  time  constant  of  the  temperature  probe  (the  upper  frequency  cutoff). 
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The  ratio  of  the  filtered  function  to  the  ideal  unfiltered  structure  function 


D(r,0,0) 


(15) 


gives  the  degradation  to  the  calculation  of  CT2  due  to  high  frequency  filtration  by 
the  probe.  Figure  2-6  shows  the  correction  required  for  different  separation 
distances  between  temperature  measurements  for  different  probe  separations. 


THERMOSONDE  PROBE  CALIBRATION 

50  msec  RESPONSE  TIME;  35  msec  AVERAGE 


1.000  2.000  3.000 

PROBE  SEPERATION  DISTANCE  (m) 

Figure  2-6:  Probe  Calibration 


3.       Computer  Program 

The  final  component  of  the  single  point  system  is  a  computer  program 
written  in  compiled  QuickBASIC  4.5  that  uses  two  converted  binary  to  ASCII  VIZ 
flight  files  as  its  input.  The  first  file  (flight. MTX)  provides  time,  temperature, 
pressure,  humidity,  spare  1,  spare  2,  spare  3,  and  spare  4  resistance  data.  The 
second  file  (flight. MDX)  provides  time  and  altitude  information. 
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A  21  point  running  mean  is  applied  twice  to  the  altitude  data  to  remove 
the  pressure  sensor  noise.  Several  different  filtering  methods  were  tried,  with  the 
double  21  point  mean  giving  the  best  results  as  shown  in  figures  2-7  through 
2-11.  Next,  the  resistance  data  from  all  four  spare  channels  were  converted  to 
temperature  by  the  method  discussed  on  page  12.  The  temperature  differences 
from  channel  1  to  the  next  data  measurement  on  channel  1  (spare  1-1)  were 
calculated  along  with  r  (the  distance  between  the  measurements)  and  Cn2  was 
computed  according  to  equation  6.  Cn2  was  next  used  to  calculate  r0  and  60, 
using  equations  8  and  9.  This  same  procedure  was  repeated  for  temperature 
differences  between  channels  2  to  1,  3  to  1,  and  4  to  1.  Summarized,  the  times 
between  temperature  measurements  where  found  as  follows: 

1)  spare  1-1  =  T1:  -  T1j+1  (=1.23  sec  apart) 

2)  spare  2-1  =  T1(  -  T2l  (=0.05  sec  apart) 

3)  spare  3-1  =  T1:  -  T3l  (=0.10  sec  apart) 

4)  spare  4-1  =  T,,  -  T4i  (=0.15  sec  apart) 
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Figure  2-8:  Effect  of  5  Point  Filter  of  Altitude  Data 
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Figure  2-9:  Effect  of  21  Point  Filter  of  Altitude  Data 
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Figure  2-10:  Effect  of  31  Point  Filter  of  Altitude  Data 
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III.  EXPERIMENTS 

A.       PROCEDURE 

1.       Maui  Flights 

The  mission  of  the  Maui  trip  was  to  aid  in  site  selection  for  a  proposed 
Air  Force,  4  meter  telescope  to  be  installed.  A  series  of  acoustic  sounder 
measurements  and  balloon  launches  occurred  at  the  Air  Force  Maui  Observation 
Station  (AMOS),  Haleakala,  Hawaii  between  26  June  to  2  July  1991  [Ref.  13]. 
Additionally,  the  1.6  meter  telescope  with  the  Itek  Compensated  Imaging  System 
(CIS)  at  AMOS  gathered  stellar  image  full  width  half  maximum  data.  The  NPS 
group  also  launched  single  point  probe  thermosonde  systems  on  28,  29,  and  30 
June  1991.  The  28  and  29  June  flights  used  the  newer  (168B)  Microsonde  II 
sondes,  while  the  30  June  flight  used  the  older  version  (Rev  X-6)  sonde. 

All  three  thermosonde  flights  were  launched  after  sunset 
(approximately  2200  local  time)  to  avoid  solar  heating  of  the  probes  [Ref.  10]. 
At  the  time  of  launch,  the  NPS  group  did  not  know  the  new  version  sondes  (28 
and  29  June)  had  different  electronic  configurations  than  the  older  version 
sondes,  therefore  the  probe  wiring  was  connected  similarly  to  the  previous 
sondes.  This  became  a  problem  because  a  resistance  voltage  divider  and  diode 
had  been  added  to  the  circuit  connected  across  spare  channel  4  to  monitor  the 
sonde  battery  voltage;  the  30  June  flight  however,  being  an  older  version  sonde 
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used  in  previous  NPS  work,  did  not  include  this  circuitry.  Although  the  resistance 
divider  and  diode  contaminated  the  results,  particularly  at  higher  altitudes  where 
the  resistances  are  in  the  100k  to  1M  range,  the  lower  altitude  data  was  not 
contaminated  excessively. 

2.       NPS  Flights 

The  NPS  group  measured  90  optically  15  November  1991  from  2130 
to  2330  local  time.  A  thermosonde  system  was  launched  from  the  roof  of  the 
Naval  Postgraduate  School's  Spanagel  Hall  during  this  same  period.  The  sonde 
was  a  new  version  Microsonde  II  with  the  voltage  divider  and  diode  removed, 
however  this  launch  included  two  probes  as  shown  in  figure  3-1. 
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Figure  3-1 :  Connecting  Probe  to  New 
Version  Sonde 


Two  thermistor  probes  connected  to  channel  1  and  channel  2  respectively  and 
separated  by  1m  provided  a  fixed  separation  distance  between  probes.  This 
provided  a  means  to  compute  horizontal  differential  temperature  measurements. 
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3.       MT.  Wilson  Flights 

In  support  of  characterization  measurements  of  Mt.  Wilson,  California, 
three  balloon  launches  occurred  in  September  1991.  None  of  these  yielded 
useful  data,  since  RF  interference  on  Mt.  Wilson  contaminated  the  measurements. 
The  Loran  and  Omega  navigational  aid  systems  used  to  determine  the  sondes' 
location  for  determination  of  wind  data  were  particularly  suspectable  to  the 
interference. 

B.       RESULTS 

1.       Maui  Flights 

Table  2  shows  the  results  of  the  thermosonde  flights  for  28,  29,  and 
30  June  1991  (one  flight  per  night)  in  comparison  with  the  acoustic  sounder  data 
and  the  AMOS  1.6  m  telescope  data  (their  measurement  of  r0  at  the  time  of  the 
balloon  flight). 

A  5  point  running  mean  of  the  temperature  data  was  used  in  the  thermosonde 
CT2  program  and  then  corrected  as  shown  on  page  13. 

TABLE  2.  CORRECTED  R0  (cm)  COMPARISON  MEASUREMENTS 


DATE 

AMOS 

NPS 

1-1 

3-1 

4-1 

AVE 

28  JUN 

10-11 

3.9 

3.1 

3.3 

3.4 

29  JUN 

7-8 

9-10 

2.26 

2.25 

2.29 

2.27 

30  JUN 

9-10 

8.9 

7.0 

7.95 
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The  flights  on  the  28  and  29  June  differed  from  the  other  systems  by  a  factor  of 
three.  This  is  believed  to  be  due  to  the  diode  circuit  in  spare  channel  4  (new 
version  sondes).  The  30  June  flight  (old  version  sonde)  shows  agreement  with 
the  AMOS  1.6  m  CIS  measurements  within  12%  discussed  in  Reference  13. 

A  problem  exists  when  computing  vertical  temperature  differences. 
Since  the  mean  atmospheric  temperature  changes  rapidly  with  altitude,  vertical 
temperature  differences  <(T(z2)  -  T(z,))2>  contain  some  unavoidable  mean  value 
changes  that  are  not  induced  by  turbulence.  This  phenomena  is  particularly 
prominent  as  the  sensors  cross  a  sharp  temperature  inversion.  This  thesis 
considered  a  variety  of  running  means  to  remove  the  T2  -  T1  mean.  Figures  3-2 
and  3-3  show  the  effect  of  5  to  21  point  running  means  on  the  computation  of  R0 
and  90. 

Figure  3-2  shows  a  monotonic  increase  in  r0  as  the  length  of  the 
running  mean  becomes  shorter.  This  demonstrates  that  there  may  be  no 
acceptable  running  mean  length  and  other  reduction  techniques  may  be 
necessary.  From  figure  3-2  it  is  apparent  that  the  channel  4-1  measurements  are 
much  less  sensitive  to  the  form  of  the  running  mean.  This  is  simply  a 
consequence  of  the  short  time  period  and  distance  between  measurements.  The 
channel  1-1  measurements  are  widely  separated  and  are  inherently  sensitive  to 
gradients  in  temperature.  In  table  2  the  channel  4-1  measurements  are  the  most 
accurate.  The  channel  1-1  measurements  in  table  2  have  been  artificially 
enhanced  by  the  5  point  running  mean  compared  to  the  21  point  results. 
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Another  source  of  systematic  errors  present  when  differencing  channel 
1  and  channels  3  or  4  are  the  fluctuations  in  time  measurements  from  channel 
1  to  3  and  1  to  4.  To  calculate  the  separation  distance  r  used  in  equation  10,  an 
average  of  100  msec  and  150  msec  time  was  used  in  the  program.  However,  the 
signal  time  occasionally  varies  by  as  much  as  10%  .  This  leads  to  noise  in  the 
separation  distance  for  channels  3-1  and  4-1  that  is  not  present  for  channel  1-1 
differencing. 
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Figure  3-2:  Running  Mean  Sensitivity 
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RUNNING  MEAN  SENSITIVITY 
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Figure  3-3:  Running  Mean  Sensitivity 

2.       NPS  Flight 

Table  3  shows  the  results  of  the  thermosonde  flight  done  on  15 
November  1991  compared  to  optical  measurements  of  the  isoplanatic  angle.  The 
dual  point  probe  and  optical  measurement  method  agreed  very  well,  while  the 
single  point  probe  are  13%  smaller  than  the  optical  measurement. 
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TABLE  3.  CORRECTED  0O  (^radians)  COMPARISON 

MEASUREMENTS 

DATE 

OPTIC 

DUAL 

1-1 

3-1 

4-1 

AVE 

15  NOV 

3-4 

3.3 

3.1 

2.2 

2.4 

2.6 

Tables  2  and  3  suggest  that  the  sonde  circuitry  may  be  inducing 
extraneous  noise  in  the  data.  The  VIZ  9000  MICROSONDE  II  sondes  measure 
resistance  by  counting  clock  pulses  of  an  oscillator  circuit  [Ref.  10].  These  clock 
pulses  (counts)  are  then  converted  to  resistance  by  the  manufacturers  software 
that  is  part  of  the  VIZ  system.  To  maintain  high  resistance  and  temperature 
resolution  (0.002  to  0.008  K),  the  sonde  alters  the  measurement  period  by  factors 
of  two  to  keep  the  counts  in  the  215  to  216  range.  It  was  found  that  this  altering 
or  "ratcheting"  happened  at  approximately  the  60  and  260  kQ  resistance 
measurements.  This  has  an  unfortunate  side  effect  on  the  small  thermistor  probe 
self  heating.  Ohmic  heating  of  the  thermistor  is  proportional  to  the  power 
dissipated  by  the  probe  as  well  as  the  time  duration  of  the  probe  current.  Figure 
3-4  shows  a  ratcheting  effect  in  the  channel  2  minus  channel  1  temperature 
differences  at  the  60  kQ  shift  point.  A  consistent  shift  of  0.5  K  occurs  above  7300 
m  in  the  temperature  differences  as  the  sonde  shifts  to  a  shorter  count  interval. 
Significant  temperature  spikes  on  the  order  of  0.05  K  occur  between  7100  and 
7300  m. 


28 


These  fluctuations  appear  to  be  the  result  of  an  instability  near  the  threshold 
point,  as  if  the  sonde  shifts  periods  back  and  forth  several  times  before 
stabilizing.  This  effect  induces  a  false  turbulence  signal  in  the  isoplanatic  angle 
on  the  order  of  5-10%. 

This  phenomena  appears  in  every  sonde  launch  investigated  to  date 
when  the  sonde  clock  circuitry  shifts  averaging  times.  Appendix  E  includes  a 
series  of  plots  that  further  demonstrates  this  effect.  Surprisingly,  the  "stutter" 
effect  appears  in  the  temperature  differences  but  not  in  the  raw  counts  recorded 
by  the  system.  It  is  clear  that  channels  5-8  each  shift  periods  at  different  times. 
The  cause  of  this  phenomenon  is  not  completely  understood  at  this  time  and 
should  be  investigated. 
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Figure  3-4:  Ratcheting  Effect  on  Delta  T 
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IV.  CONCLUSIONS  AND  RECOMMENDATIONS 

A.       CONCLUSIONS 

The  new  high  speed  5  mil  VECO  thermistor  probe  proved  to  be  a  large 
improvement  in  the  single  point  probe  method  of  measuring  atmospheric 
turbulence.  The  results  of  measurements  of  r0  and  90  are  now  within  10  to  12 
percent  of  the  other  methods  of  measuring  the  turbulence. 

The  new  version  sondes  are  more  useful  than  the  older  version  by  allowing 
the  use  of  all  4  spare  channels  for  data  collection.  The  noise  in  the  new  version 
sonde  is  consistent  with  what  was  measured  in  the  old  version,  therefore  the 
noise  filtering  methods  used  for  the  old  sondes  are  the  same  for  the  new  version 
sondes.  A  successful  method  of  filtering  the  altitude  data  was  found  to  be  a 
double  21  point  running  mean. 

This  thesis  found  that  a  running  mean  applied  to  first  order  temperature 
differences  was  not  sufficient  for  trend  removal.  Coherence  length  computed 
with  5-21  point  running  means  showed  a  monotonic  increase  for  shorter  running 
means.  It  appears  that  no  running  mean  is  suitable.  Other  data  reduction 
techniques  will  be  necessary. 
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Systematic  spikes  in  Microsonde  II  temperature  data  were  found  near  the 
resistance  range  change  points.  These  temperature  spikes  contaminated  the 
isoplanatic  angle  by  10%.  They  will  have  to  be  removed  by  software  or  hardware 
modifications  in  future  work. 

The  improved  high  speed  thermistor,  along  with  the  new  version  sonde  and 
a  computer  program  incorporating  altitude  filtering  and  high  order  temperature 
differencing  will  allow  the  NPS  group  to  have  a  working  single  point  probe 
system.  This  system  has  the  ability  to  plot  altitude  verses  Cn2  displaying  where 
the  layers  of  optical  turbulence  are  encountered. 


B.       RECOMMENDATIONS 

To  gain  more  confidence  in  the  filtering  methods  used  in  the  computer 
program  and  to  establish  a  statistical  data  base  for  measurement  comparisons, 
a  number  of  flights  should  be  launched  with  the  two  point  probe  system.  These 
obviously  should  cover  a  wide  range  of  both  "good  seeing"  and  "poor  seeing" 
days.  Another  advantage  of  doing  numerous  flights  to  establish  a  data  base  is 
that  a  calibration  correction  factor  could  then  be  applied  to  the  3-1  and  4-1 
channels. 
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A  continued  study  of  how  to  handle  the  temperature  inversions  with  a  signal 
filter  needs  to  be  carried  out  and  incorporated  into  the  computer  program.  This 
is  probably  the  future  research  area  that  would  produce  the  greatest 
improvement  to  the  system. 

It  is  unknown  if  the  variations  in  the  timing  intervals  between  channel  1  to 
4  can  be  handled  or  removed  well  enough  to  bring  into  agreement 
measurements  across  all  channels.  If  the  noise  is  not  improved,  then  output 
measurements  should  be  restricted  to  differencing  channel  1-1  for  a  single  point 
probe  or  2-1  for  a  dual  point  probe. 

Shortening  the  time  between  scans  is  a  hardware  improvement  to  the  VIZ 
Microsonde  II  that  will  improve  performance.  If  this  scan  interval  could  be 
shortened  from  the  present  1 .23  sec  to  approximately  0.5  sec  with  the  same  low 
noise,  it  would  be  a  big  improvement  for  use  as  a  single  point  system.  This  is 
because  0.15  sec  tends  to  be  too  short  for  a  50  ms  probe  response  time,  while 
1.23  sec  is  too  long  (which  corresponds  to  =5  m  spacing).  Therefore  0.5  sec 
(~2  m  spacing)  should  improve  the  measurements. 
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APPENDIX  A.  PROGRAM  FOR  COMPUTING  Cnz,  R0,  AND  0O 


DECLARE  SUB  Crout  (a!(),  b!(),  c!) 
DECLARE  FUNCTION  temp!  (r!,  b!()) 
DEFINT  I-N 
OPTION  BASE  1 
DIM  xm(0  TO  10) 
DIM  a(4,  4) ,  b(4) 

This  program  reads  the  ASCII   .mtx  and  .mdx  file  created  by 

the  VIZ  9000  SYSTEM  and  performs  a  least  squares  fit 

creating  a  temperature  data  output. 

It  then  filters  the  altitude  data  and  computes  CN~2,  Ro,  and 

theta_o  for  4  spare  channels  of  data  from  the  temperature 

data  creating  an  altitude  verses  CN^2  '  output  file  to  be 

plotted  with  SUPERCALC. 

DLW  Sep  11  1990   (temperature  conversion  program) 

CRH  Nov  10  1991 

To  declare  dynamic  data  array,  invoke  Quickbasic  with  /ah 

option 


LIST 
start 

stp  = 


OF  VARIABLES  AND  CONSTANTS: 
=  starting  altitude  data  point 
stopping  altitude  data  point 


dta 
dta 

dta 
dta 

dta 
dta 
dta 
dta 
dta 
dta 
dta 
dta 
dta 
dta 
dta 
dta 
dta 
b(l 

b(2) 


, 1) =  time  (sec) 

,2)=  pressure  (mbar) 

, 3)=  VIZ  temp  (C) 

,4)=  humidity  (%) 


, 5) =  spare 

,6)=  spare 

,7)  =  spare 

,8)=  spare 


rol=  Ro  for  spare  1-1 
Ool=  theta_o  for 

spare  1-1 
ro2  =  Ro  for  spare  2-1 
Oo2=  theta_o  for 
spare  2-1 
ro3=  Ro  for  spare  3-1 
Oo3=  theta_o  for 
spare  3-1 
temp  (K) 
temp  (K) 
temp  (K) 
temp  (K) 


resistance 

resistance 

resistance 

resistance 
,9)=  VIZ  altitude  (m) 
,10)=  first  21  pt  mean  altitude 
,11)=  2nd  21  pt  mean  altitude 
,12)=  delta  T  for  all  spares  (reused) 
,13)=  mean  delta  T  for  all  spares  (reused) 
,14)=  CN"2  for  spare  1-1  ro4=  Ro  for  spare  4-1 

,15)=  CN~2  for  spare  2-1         0o4=  theta_o  for 
,16)=  CN~2  for  spare  3-1  spare  4-1 

,17)=  CN*2  for  spare  4-1 
coefficient  of  polynomial  to  convert  resistance  to 
temperature 
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'  b(3)=  coefficient  of  polynomial  to  convert  resistance  to 

temperature 
'  b(4)=   ,, 

1  low  =  number  of  data  pts  lost  due  to  running  mean  size 
'  high  =  number  of  data  pts  lost  due  to  running  mean  size  at 

end  of  flight  for  spare  1-1  only. 
•  ST  =  desired  starting  altitude 
'  RM  =  desired  running  mean  size 
1  pts  =  number  of  vertical  data  points 


REM  $DYNAMIC 

DIM  dta(6000,  17) 

tzero  =  273.15  '  zero  Celcius 

INPUT  "enter  MTX  filename  ";  MTX$ 

PRINT 

PRINT  MTX$ 

OPEN  MTX$  FOR  INPUT  AS  #1 

INPUT  "enter  MDX  filename  ";  MDX$ 

PRINT 

PRINT  MDX$ 

OPEN  MDX$  FOR  INPUT  AS  #2 

1  =  0 

DO  UNTIL  E0F(1) 

1  =  1  +  1 

INPUT  #1,  dta(l,  1),  dta(l,  2),  dta ( 1 ,  3),  dta ( 1 ,  4), 
dta(l,  5),  dta(l,  6),  dta ( 1 ,  7),  dta ( 1 ,  8) 
LOOP 

1  =  0 

DO  UNTIL  EOF (2) 

1  =  1  +  1 

INPUT  #2,  x,  y,  dta(l,  9)   'x  and  y  are  unneeded  variables 
LOOP 

'FIND  START  AND  STOP  OF  FLIGHT 

N  =  1 

FOR  i  =  1  TO  (N  /  2) 

IF  (dta(i,  9)  <  dta ( i  +  1,  9))  AND 

(dta(i  +  1,  9)  <  dta(i  +  2,  9))  AND 

(dta(i  +  2,  9)  <  dta(i  +  3,  9))  AND 

(dta(i  +3,  9)  <  dta(i  +  4,  9))  AND 

(dta(i  +  1,  9)  -  dta(i,  9))  >  2  AND 
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(dta(i  +  2,  9)  -  dta(i  +  1,  9))  >  2 
THEN  GOTO  10 
NEXT  i 

10  start  =  i      'exited  loop  when  conditions  were  met 

FOR  i  =  N  TO  (N  /  2)  STEP  -1 

IF  dta(i,  9)  >  dta ( i  -  1,  9)  AND 

dta(i  -  1,  9)  >  dta(i  -  2,  9)  AND 
dta(i  -  2,  9)  >  dta(i  -  3,  9)  THEN  GOTO  2  0 
NEXT  i 
20  stp  =  i        'exited  loop  when  conditions  were  met 

CLOSE  #1:  CLOSE  #2 

•CORRECT  ALTITUDE  (VIZ  MAX  ALTITUDE  INPUT  IS  2500  m) 

INPUT  "enter  altitude  correction  (m)";  zcor 

PRINT 

PRINT  zcor 

FOR  i  =  start  TO  stp 

dta(i,  9)  =  dta(i,  9)  +  zcor 
NEXT  i 

'DO  DOUBLE  21  PT  RUNNING  MEAN  ON  ALTITUDE  (Z) 

'first  21  point  mean 
FOR  i  =  start  TO  (stp  -  20) 
num!  =  0! 
FOR  j  =  i  TO  (i  +  20) 

num!  =  num!  +  dta(j,  9) 
NEXT  j 
dta(i  +  10,  10)  =  num!  /  21! 
NEXT  i 

'second  21  point  mean 
FOR  i  =  (start  +  10)  TO  (stp  -  30) 
num !  =  0 ! 
FOR  j  =  i  TO  (i  +  20) 

num!  =  num!  +  dta(j,  10) 
NEXT  j 
dta(i  +  10,  11)  =  num!  /  21! 
NEXT  i 

'  compute  average  ascent  distance  between  data  points  at  start 
'  of  flight 

zlave!  =  (dta(start  +  24,  11)  -  dta(start  +  20,  11))  /  4! 
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'  compute  average  ascent  distance  between  data  points  at  end 
1  of  flight 

z2ave!  =  (dta(stp  -  20,  11)  -  dta(stp  -  24,  11))  /  4! 

'  compute  averaged  altitudes  to  fill  in  first  20  data  pts  lost 

•  by  double  filter 

FOR  j  =  1  TO  2  0 

dta(j,  11)  =  dta(start  +  20,  11)  -  ((21  -  j)  *  (zlave!)) 
NEXT  j 

j  =  1 

FOR  i  =  start  TO  (start  +  19) 

dta(i,  11)  =  dta(j,  11) 

j  =  j  +  1 
NEXT  i 

'  compute  averaged  altitudes  to  fill  in  last  20  data  pts  lost 

•  by  double  filter 

FOR  j  =  1  TO  2  0 

dta(j,  11)  =  dta(stp  -  20,  11)  +  (j  *  z2ave!) 
NEXT  j 

j  =  1 

FOR  i  =  (stp  -  19)  TO  stp 

dta(i,  11)  =  dta( j ,  11) 

j  =  j  +  1 
NEXT  i 

'  The  following  is  a  program  to  convert  resistance  data  to 

•  temp  data 

1  begin  least  sguares  fit 

m  =  4 

FOR  i  =  1  TO  m 

FOR  j  =  1  TO  m 
a(i,  j)  -  0! 

NEXT  j 

b(i)  =  0! 
NEXT  i 

m2  =  m  +  m  -  2 

FOR  k  =  start  TO  stp 

xk  =  LOG(dta(k,  5) ) 

yk  =  1!  /  (dta(k,  3)  +  tzero) 

xm(0)  =  1! 


37 


FOR  i  =  1  TO  m2 

xm(i)  =  xk  *  xm(i  -  1)  '  xm  =  xAm 

NEXT  i 

FOR  i  =  1  TO  m 
±2    =    i  -  2 

FOR  j  =  i  TO  m 

a(i,  j)  =  a(i,  j)  +  xm(i2  +  j)  '  x~(i-l+j-l) 

NEXT  j 

b(i)  =  b(i)  +  yk  *  xm(i  -  1) 
NEXT  i 

NEXT  k 

'  set  lower  half  of  matrix 

FOR  i  =  1  TO  m 

FOR  j  =  1  TO  i 

a(i,  j)  =  a(j,  i) 

NEXT  j 
NEXT  i 

c  =  m 

CALL  Crout(a() ,  b() ,  c) 

PRINT  b(l),  b(2)  ,  b(3),  b(4) 

1  convert  to  temperature 

FOR  i  =  start  TO  stp 

dta(i,  3)  =  dta(i,  3)  +  tzero 

resistance  =  dta(i,  5) 

dta ( i ,  5)  =  temp (resistance,  b()) 

resistance  =  dta ( i ,  6) 

dta ( i ,  6)  =  temp (resistance,  b()) 

resistance  =  dta(i,  7) 

dta(i,  7)  =  temp (resistance,  b()) 

resistance  =  dta(i,  8) 

dta(i,  8)  =  temp (resistance,  b()) 
NEXT  i 

•RENUMBER  ARRAY  FROM  1  TO  N 

FOR  i  =  start  TO  stp 

FOR  j  =  1  TO  11 

dta(i  +  1  -  start,  j)  =  dta(i,  j) 

NEXT  j 
NEXT  i 

N  =  (i  -  start  +  1) 
PRINT 
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PRINT  "starting  altitude  =  ";  dta ( 1 ,  11) 
PRINT  "stopping  altitude  =  ";  dta(N,  11) 

1  allow  user  to  select  new  starting  altitude  to  examine 
1  portion  of  flight 

PRINT 

INPUT  "enter  new  desired  starting  altitude";  ST 

PRINT 

PRINT  ST 

FOR  i  =  1  TO  N 

IF  dta(i,  11)  >  ST  THEN  GOTO  30 
NEXT  i 
30  start  =  i 
PRINT 
PRINT  "New  starting  altitude  =  ";  dta (start,  11) 

•  compute  number  of  data  points 

stp  =  N 

pts  =  stp  -  start  +  1 

PRINT 

PRINT  "number  of  data  points  are  ";  pts 

'ENTER  RUNNING  MEAN  TO  BE  USED  FOR  dT 

INPUT  "enter  size  of  running  mean  (odd  number)  for  temp";  RM 

PRINT 

PRINT  "dT  running  mean  is  ";  RM 

•COMPUTE  CT~2  FOR  SPARE  1  -  SPARE  1  TEMP  DIFFERENCE 

'COMPUTE  DELTA  T 

FOR  i  =  start  TO  (stp  -  1) 

dta(i,  12)  =  dta(i,  5)  -  dta ( i  +  1,  5) 
NEXT  i 

'DO  RUNNING  MEAN  OF  dT 

low  =  (RM  -  1)  /  2 

high  =  (RM  +  1)  /  2 

FOR  i  -  start  TO  (stp  -  RM) 

num!  =  0! 

FOR  j=iTO(i+RM-l) 

num!  =  num!  +  dta(j,  12) 

NEXT  j 

dta(i  +  low,  13)  =  num!  /  RM! 
NEXT  i 
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'compute  cn2 

FOR  i  =  (start  +  low)  TO  (stp  -  high) 
dta(i,  14)  =  (((.000079  *  dta  ( i ,  2))  /  (dta(i,  3)  A  2))  A  2)* 
((dta(i,  12)  -  dta(i,  13))  A  2)  /  ((dta(i  +  1,  11) 
-  dta(i,  11) )  A  ( .6667) ) 
NEXT  i 

h!  =  2.1  *  (1.46  *  (2  *  3.14159  /  500E-09)  A  2)  A  -0.6 
sum!  =  0! 

'  compute  Ro  (CN~2  *  dz) 

FOR  i  =  (start  +  low)  TO  (stp  -  high) 
sum!  =  sum!  +  ((dta(i  +  1,  11)  -  dta(i,  11))  *  dta ( i ,  14)) 

•dz=dta(i+l, 11) -dta (i, 11) 
NEXT  i 

sum2 !  =  sum!  A  (- . 6) 
rol!  =  h!  *  sum2 ! 

'COMPUTE  Oo  for  1-1  (CN~2  *  dz  *  zA5/3) 

f!  =  (2.905  *  (2  *  3.14159  /  500E-09)  A  2)  A  -0.6 

sumO!  =  0! 

FOR  i  =  (start  +  low)  TO  (stp  -  high) 

sumO!  =  sumO!  +  (((dta(i  +  1,  11)  -  dta(i,  11))  * 
dta(i,  14))  *  (dta(i,  11)  -  dta(start,  11))  A  1.6667) 
NEXT  i 

sum02 !  =  sumO!  A  (-.6) 
Ool!  =  f!  *  sum02 ! 

'COMPUTE  CT~2  FOR  SPARE  2  -  SPARE  1  TEMP  DIFFERENCE 

'compute  delta  T 

FOR  i  =  start  TO  stp 

dta(i,  12)  =  dta(i,  5)  -  dta(i,  6) 
NEXT  i 

'do  running  mean  of  dT 

FOR  i  =  start  TO  (stp  -  RM  +  1) 

num !  =  0 ! 

FOR  j=iTO(i+RM-l) 

num!  =  num!  +  dta(j,  12) 

NEXT  j 

dta(i  +  low,  13)  =  num!  /  RM! 
NEXT  i 

' compute  cn2 
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FOR  i  =  (start  +  low)  TO  (stp  -  low) 

dta(i,  15)  =  (((.000079  *  dta(i,  2))  /  (dta(i,  3)  A  2))  "2)* 
((dta(i,  12)  -  dta(i,  13))  A  2)  /  ((((dta(i  +  1,  11) 
-  dta(i,  11))  /  1.23)  *  .05)  A  (.6667)) 
NEXT  i 
sum!  =  0! 
sum2 !  =  0! 
FOR  i  =  (start  +  low)  TO  (stp  -  low) 

sum!  =  sum!  +  (dta(i,  15)  *  (dta(i  +  1,  11)  -  dta(i,  11))) 
NEXT  i 

sum2 !  =  sum!  *  (-.6) 
ro2 !  =  h!  *  sum2 ! 

'compute  Oo  for  2-1 

sumO !  =  0 ! 

FOR  i  =  (start  +  low)  TO  (stp  -  low) 

sumO!  =  sumO!  +  (((dta(i  +  1,  11)  -  dta ( i ,  11))  * 

dta(i,  15))  *  (dta(i,  11)  -  dta(start,  11))  ~  1.6667) 
NEXT  i 

sum02 !  =  sumO!  "     (-.6) 
Oo2!  =  f!  *  sum02! 

'COMPUTE  CT~2  FOR  SPARE  3  -  SPARE  1  TEMP  DIFFERENCE 

'compute  dT 

FOR  i  =  start  TO  stp 

dta(i,  12)  =  dta(i,  5)  -  dta(i,  7) 
NEXT  i 

•compute  running  mean  of  dT 

FOR  i  =  start  TO  (stp  -  RM  +  1) 
num!  =  0! 
FOR  j  =  i  TO  (i  +  RM  -  1) 

num!  =  num!  +  dta(j,  12) 
NEXT  j 
dta(i  +  low,  13)  =  num!  /  RM! 
NEXT  i 

'compute  cn2 

FOR  i  =  (start  +  low)  TO  (stp  -  low) 

dta(i,  16)  =  (((.000079  *  dta  ( i ,  2))  /  (dta(i,  3)  A  2))  "2)* 
((dta(i,  12)  -  dta(i,  13))  A  2)  /  ((((dta(i  +  1,  11) 
-  dta(i,  11))  /  1.23)  *  .1)  ~  (.6667))     » dz / 1 . 2 3=  vel ' y 

' r=vel ' y  *  0 . lsec 
NEXT  i 
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sum!  =  0! 

sum2  !  =  0! 

FOR  i  =  (start  +  low)  TO  (stp  -  low) 

sum!  =  sum!  +  (dta(i,  16)  *  (dta(i  +  1,  11)  -  dta  ( i ,  11))) 
NEXT  i 

sum2 !  =  sum!  *  (-.6) 
ro3 !  =  h!  *  sum2 ! 

'compute  Oo  for  3-1 

sumO!  =  0! 

FOR  i  =  (start  +  low)  TO  (stp  -  low) 

sumO!  =  sumO!  +  (((dta(i  +  1,  11)  -  dta(i,  11))  * 

dta(i,  16))  *  (dta(i,  11)  -  dta (start,  11))  A  1.6667) 
NEXT  i 

sum02 !  =  sumO!  A  (-.6) 
Oo3!  =  f!  *  sum02 ! 

'COMPUTE  CT~2  FOR  SPARE  4  -  SPARE  1  TEMP  DIFFERENCE 

'compute  delta  T 

FOR  i  =  start  TO  stp 

dta(i,  12)  =  dta(i,  5)  -  dta(i,  8) 
NEXT  i 
FOR  i  =  start  TO  (stp  -  RM  +  1) 

num!  =  0! 

FOR  j  =  i  TO  (i  +  RM  -  1) 

num!  =  num!  +  dta(j,  12) 
NEXT  j 

dta(i  +  low,  13)  =  num!  /  RM! 
NEXT  i 

'compute  cn2 

FOR  i  =  (start  +  low)  TO  (stp  -  low) 

dta(i,  17)  =  (((.000079  *  dta  ( i ,  2))  /  (dta(i,  3)  *  2))  "2)* 
((dta(i,  12)  -  dta(i,  13))  A  2)  /  ((((dta(i  +  1,  11) 
-  dta(i,  11))  /  1.23)  *  .15)  A  (.6667)) 
NEXT  i 
sum!  =  0! 
sum2 !  =  0 ! 
FOR  i  =  (start  +  low)  TO  (stp  -  low) 

sum!  =  sum!  +  dta ( i ,  17)  *  (dta(i  +  1,  11)  -  dta ( i ,  11)) 
NEXT  i 

sum2 !  =  sum!  A  (-.6) 
ro4 !  =  h!  *  sum2 ! 

'COMPUTE  Oo  for  4-1 
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sumO!  =  0! 

FOR  i  =  (start  +  low)  TO  (stp  -  low) 

sumO!  =  sumO!  +  (((dta(i  +  1,  11)  -  dta  ( i ,  11))  * 

dta(i,  17))  *  (dta(i,  11)  -  dta(start,  11))  *    1.6667) 
NEXT  i 

sum02 !  =  sumO!  A  (-.6) 
Oo4 !  =  f!  *  sum02! 

•OUTPUT  TO  A  FILE  ALTITUDE  VERSES  CN2  FOR  EACH  CHANNEL  TO  BE 
'PLOTTED  WITH  SUPERCALC 

INPUT  "enter  output  file  name  ",  output$ 
OPEN  output$  FOR  OUTPUT  AS  #3 

FOR  i  =  start  TO  stp 

WRITE  #3,  dta(i,  11),  dta ( i ,  14),  dta ( i  ,  15),  dta ( i ,  16), 
dta(i,  17) 
NEXT  i 


PRINT 

PRINT 

"Roll  = 

-  ii 

,  rol! 

PRINT 

"Ro2  = 

ro2  ! 

PRINT 

"Ro3  = 

ii 

ro3  ! 

PRINT 

"Ro4  = 

ii 

ro4  ! 

PRINT 

"Ool  = 

ii  . 

Ool! 

PRINT 

"Oo2  = 

ii  . 

i 

Oo2! 

PRINT 

"Oo3  = 

ii  . 

r 

Oo3  ! 

PRINT 

"Oo4  = 

ii  . 

/ 

Oo4  ! 

END 


REM  $STATIC 

DEFSNG  I-N 

SUB  Crout  (a() ,  b() ,  c) 
DEFINT  I-M 

'  this  subroutine  solves  the  set  of  simultaneous  equations 
'  A*X=B  where  A  is  a  symmetric  matrix  using  the  Crout  or 


'  Cholesky  procedure. 

'  The  original  A  and  B  arrays  are  destroyed 

m  =  c 

all  =  a(l,  1) 

IF  all  =  0  THEN  GOTO  singular 

FOR  i  =  2  TO  m 

a(l,  i)  =  a(i,  1)  /  all 
NEXT  i 
b(l)  =  b(l)  /  all 
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'set  up  auxiliary  matrix 

FOR  j  =  2  TO  m 
jl  =  j  -  1 
FOR  i  =  j  TO  m 

asave  =  0! 

FOR  k  =  1  TO  jl 

asave  =  asave  +  a(i,  k)  *  a(k,  j) 

NEXT  k 

a(i,  j)  =a(i,  j)  -  asave 

IF  i  >  j  THEN  a(j,  i)  =  a(i,  j)  /  a(j,  j) 

NEXT  i 

bsav  =  0! 

FOR  k  =  1  TO  jl 

bsav  =  bsav  +  a(j,  k)  *  b(k) 

NEXT  k 

ajj  =  a(j,  j) 

IF  ajj  =  0  THEN  GOTO  singular 

b(j)  =  (b(j)  -  bsav)  /  ajj 

NEXT  j 
i 

'  back  substitution  and  calculate  determinant 
i 

ml  =  m  -  1 
FOR  i  =  1  TO  ml 
bsav  =  0! 
mi  =  m  -  i 
mil  =  mi  +  1 
FOR  j  =  mil  TO  m 

bsav  =  bsav  +  a(mi,  j)  *  b(j) 
NEXT  j 

b(mi)  =  b(mi)  -  bsav 
NEXT  i 
GOTO  finished 

singular : 

PRINT  "***Matrix  is  singular***" 
finished : 
END  SUB 

DEFSNG  I-M 

FUNCTION  temp  (r,  b()) 

lnr  =  LOG(r) 

temp  =  1!  /  (b(l)  +  lnr  *  (b(2)  +  lnr  *  (b(3)  +  b(4)  *  lnr))) 

END  FUNCTION 
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APPENDIX  B.  NOTES  ON  OPERATION  OF  VIZ-9000  SYSTEM 

The  following  are  notes  to  support  the  operation  of  the  VIZ  system.    They 
are  to  be  used  in  conjunction  with  the  manufacturer's  manual.    These  notes  are 
a  brief  summary  of  the  steps  the  NPS  group  found  to  be  important  upon  initial 
indoctrination  of  how  to  use  the  VIZ  system. 
PART  I 
•  The  initial  wiring  of  the  system  should  be  as  shown  in  figure  B-1. 


COMPUTER 


ANT       403  REC       0M      L0R      TK1    J   0M       L0R     TK2      LNCH     INTF 
REM     REM  /  LOC     LOC 


Figure  B-1:  Wiring  Diagram 
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PART  II 

•  After  properly  connecting  all  the  cables,  and   properly  installing  the 

navigation  and  receiver  antennas,  the  system  should  be  booted  up  ensuring  the 

"rack"  is  turned  on  before  the  computer. 

(WARNING:  be  sure  the  nav  antenna  is  not  too  close  to  the 
receiver  antenna  before  starting  the  system.    The 
receiver  antenna  makes  a  complete  rotation  upon 
energizing  and  can  damage  a  nearby  nav  antenna.) 

•  The  next  step  is  to  run  "SETUP".     SETUP  consists  of  8  parts.     The 

following  configuration  should  be  used  for  the  Monterey  Bay  area: 

***Whenever  data  is  selected  or  entered,  to  leave  the  present  screen 
and  return  to  the  previous  menu,  press  "escape"  and  follow  the 
directions  on  the  screen. 

1)  Station  Installation: 

Customer  Name:    NPS 

Station  ID:    NPS 

LAT:  36°  35'  42.0  N 

LONG:  121°  52'  57.0  W 

WMO  station  no:  0 

Station  Ht:  38  m    (2500  m  MAX) 

Local  time:    0 

Surface  data  timer:    0 

Sonde  type:  2 

PTU  selection:    Y 

2)  Receiver  Installation: 

start-stop:  395-410  hz 

3)  Navaid  Installation: 

Zeemet  system:  LORAN 
Tracker  1 :  LORAN 
Tracker  2:  OMEGA 
****lf  desire  an  omega  launch,  then  just  toggle  the  zeemet  system 
to  omega  (by  using  space  bar),   nothing  else  has  to  be 
changed. 


46 


****For  cross  chain  loran  flight,  change  tracker  2  to  loran. 

After  entering  "LORAN,  LORAN,  OMEGA",  go  to  tracker  1  and  press  "enter".   This 

will  open  a  sub-window  to  choose  what  loran  chain  you  want.    Highlight  "9940" 

and  press  the  space  bar  so  the  "check  mark"  moves  to  9940.    Next  press 

"escape",  go  to  tracker  2,  press  "enter",  and  now  select  for  omega: 

OMEGA 
Japan:  include 

Norway:  exclude 

Liberia:  exclude 

Hawaii:  include 

North  Dakota:     common 
La  Reunion:        include 
Argentina:         include 
Australia:         include 

4)  Wind  Computation  (done  in  height) 

LORAN 

1.0  60  2  15 

2.  1000      90  3  15 

3.  10000     120  4  15 

OMEGA 

1.  0           120         4  15 

2.  1000       180         5  15 

3.  10000     240         7  15 

5)  Loran  Chain  File: 

-leave  untouched 

6)  Antenna  Configuration: 

The  cycle  time  needs  to  be  rechecked  at  least  yearly.    Presently  the  cycle  time 

is  48.3  sec. 

****Remember  at  the  end  of  this  menu  to  adjust  antenna  center  for  each 
site  (need  magnetic  variation  for  site).    For  Monterey,  magnetic 
variation  is  16°  East,  therefore  compass  heading  should  read  074 
(so  yagi  antenna  points  090  T). 
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6)  Antenna  Configuration  (continued): 

mode:  auto 

select:  omni 

systems:  both 

switch:  4000  m  down  range 

When  this  information  is  entered,  antenna  will  cycle  and  stop  at  its  midpoint. 

This  is  when  the  orientation  should  be  adjusted  (adjust  antenna  to  most  probable 

wind  azimuth). 

PART  III 

•  To  run  a  flight: 

1)  choose  which  type  of  sonde  (loran  or  omega) 

2)  fill  balloon  and  attach  sonde 

3)  get  permission  from  Monterey  Airport  tower  to  launch  (375-1214) 

4)  type  "START" 

•  When  menu  appears,  enter  sonde  serial  number  and  flight  ID: 

flight  ID  =  month  +  day  +  hour  (2  digits  each) 
EX:  031121   =  March  11,  2100 

•  Once  this  is  entered,  check  Navaid  status.  Use  instruments  to  enter  surface 
data. 

•  Once  sonde  is  prepared,  check  receiver  status  (signal  strength  should  be  200- 
255).  Receiver  should  lock  on  sonde  after  a  few  minutes.  If  not,  go  to  receiver 
display  [Alt-s]  to  observe  sweep  of  spectrum. 

•  Next  check  calibration.  Check  to  see  that  temperature  and  humidity  lock- 
in  numbers  match  those  on  sonde  sticker  (if  close,  go  with  computer,  if  way  off, 
go  with  sonde  label).    Enter  offsets  if  needed. 


48 


•  Finally,  when  "?'s"  on  1-4  have  changed  to  "*'s",  arm  for  launch.  You  can 
launch  the  sonde  at  anytime  after  you  press  arm  for  launch.  However,  do  not 
hold  the  sonde  too  long  because  you  only  have  two  hours  of  allotted  time  from 
the  "arm  for  launch"  time  for  a  complete  flight. 

•  Wait  a  few  minutes  for  PTU  data  to  show  up  on  the  screen. 

•  Once  you  have  PTU  data,  go  to  "Data  Display  Mode".  Choose  the  PTU 
plot  display. 

•  Check  ascent  rate  for  180-240  m/min. 

•  At  any  time  in  the  flight  you  can  check  on  the  pointing  of  the  antenna  and 

S/N  ratio  levels  by  going  to  "utilities". 

****Do  not  go  to  receiver  mode  in  utilities  once  launch  has  started.    This 
will  cause  a  stopping  of  data  collection. 

•  When  flight  has  finished,  you  should  be  in  either  PTU  plot  or  general  flight 
menu  screen.    Now  go  to  utilities  to  save  flight  data  to  disk. 

•  After  this  is  done,  you  can  exit  flight  menu. 
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APPENDIX  C.  PLOTS  OF  C, 
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Figure  C-1:  Altitude  verses  Cn2 
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FLIGHT  AF062922 
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Figure  C-2:  Altitude  verses  Cn2 
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FLIGHT  AF0G:i021 
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Figure  C-3:  Altitude  verses  Cn2 
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Figure  C-4:  Altitude  verses  Cn2 
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FLIGHT  NS111522 
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Figure  C-5:  Altitude  verses  Cn5 
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APPENDIX  D.  PROBE  HEATING 

Figure  D-1  is  a  plot  showing  the  effect  of  probe  heating  from  channel  1  to 
4.  As  shown,  the  probe  temperature  increases  from  channel  1  to  2,  then 
decreases  on  channel  4.  This  is  because  the  probe  is  off  between  channels  2 
and  4  (no  data  on  channel  3)  allowing  the  probe  to  cool. 


FLIGHT  AF063021 

ALTITUDE  VS  PROBE  TEMPERATURE 

3,460 

3,450 

"g  3,440 

UJ 

§3,430 

h 

<  3,420 

3,410 

3,400 
2c 

♦ 

*          :  + 

*           + 

♦ 

+ 
*            + 

»           .+ 

»             + 
•             + 
*           + 

CHANNEL  NO. 

•  CHANNEL  1 
+  CHANNEL  2 

*  CHANNEL  4 

12           282.1          282.2         282.3         28 
PROBE  TEMPERATURE  (K) 

2.4 

Figure  D-1.  Probe  Heating 
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APPENDIX  E.  EFFECTS  OF  CONVERTING  COUNTS  TO  RESISTANCE 

The  VIZ  Microsonde  II  measures  resistance  by  counting  clock  pulses  and 
then  converting  these  counts  into  resistance  values.  In  order  to  keep  resolution 
high,  the  counts  stay  between  214  and  216.  When  the  counts  fall  below  214  the 
system  doubles  the  number  of  periods  counted  (racheting  effect).  Figures  E-1 
through  E-7  show  the  effects  of  these  rachet  points.  The  two  times  the  system 
racheted  where  found  to  be  at  resistances  of  60  and  260  kQ. 
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Figure  E-1 :  Racheting  Effect  on  Counts 
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Figure  E-2:  Racheting  Effect  on  Counts 
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Figure  E-3:  Racheting  Effect  on  Counts 
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Figure  E-4:  Racheting  Effect  on  Counts 
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Figure  E-5:  Racheting  Effect  on  Delta  T 


60 


FLIGHT  AF063021 

ALTITUDE  VS  DELTA  T 
CHANNEL  4-1 


r 
C 
\s 

U 

D 

h 

j 
< 


■ 


DELIA 


Figure  E-6:  Racheting  Effect  on  Delta  T 
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Figure  E-7:  Racheting  Effect  on  Processed  Data 
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